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Abstract Wc discuss the diagnostic potential of high cadence ultraviolet spec- 
tral data when transient ionization is considered. For this we use high cadence 
UV spectra taken during the impulsive phase of a solar flares (observed with 
instruments on-board the Solar Maximum Mission) which showed excellent cor- 
respondence with hard X-ray pulses. The ionization fraction of the transition 
region ion O v and in particular the contribution function for the O v 1371 A 
line are computed within the Atomic Data and Analysis Structure, which is a 
collection of fundamental and derived atomic data and codes which manipulate 
them. Due to transient ionization, the O v 1371 A line is enhanced in the first 
fraction of a second with the peak in the line contribution function occurring 
initially at a higher electron temperature than in ionization equilibrium. The 
^ , rise time and enhancement factor depend mostly on the electron density. The 

fractional increase in the O v 1371 A cmissivity due to transient ionization 
can reach a factor of 2-4 and can explain the fast response in the line flux of 
transition regions ions during the impulsive phase of flares solely as a result 
' of transient ionization. This technique can be used to diagnostic the electron 

■ temperature and density of solar flares observed with the forth-coming Interface 

, Region Imaging Spectrograph 
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1. Introduction 



The solar atmosphere contains a range of highly dynamic features, e.g. flares, 
jets, spicules, etc. In order to fully observe these events, requires high cadence, 
high spatial and high spectral resolution instruments. However, due to vari- 
ous scientific and technical considerations, trade-offs in the instrument design 
must be made by the investigators. Milli-sccond data is commonly acquired in 
the radio and hard X-rays frequencies, while at ultraviolet wavelengths, high 
cadence is generally not available with most ultraviolet spectrographs. The 
highest cadence UV spectrometer to date was the Ultraviolet Spectrometer and 
Polarimeter (UVSP) which was flown on the Solar Maximum Mission (SMM) 
in the early 1980's. Several flares were observed simultaneously in the UV at 
selected individual spectral lines, e.g. O v 1371 A, Si iv 1394 A etc, and at hard 
X-ray energies of several keV. Several papers reported simultaneous (to within 
a fraction of a second) increases in the UV lines and at hard X-ray energies 
(Cheng et al. 1981, Woodgate et al. 1983, Poland et al. 1984, Cheng et al. 1984). 
A more in-depth study by Cheng et al. (1988) found that bursts observed in the 
O v 1371 A line lagged the hard X-ray bursts by only 0.3 s to 0.7 s. No firm 
explanation was given for such a fast rise in the UV line flux. 

None of the various spectrographs launched since SMM has had the ability to 
obtain high cadence spectral data. However, that will shortly change with the 
launch of the Interface Region Imaging Spectrograph (IRIS j 1 ] which is centered 
on four themes relating to solar physics and space weather, namely: Which types 
of non-thermal energy dominate in the chromosphere and beyond, how does the 
chromosphere regulate mass and energy supply to the corona and hcliosphcrc and 
how does magnetic flux and matter rise through the lower atmosphere, and what 
role docs flux emergence play in flares and mass ejections? The spectrograph will 
have the ability to obtain spectra with a cadence down to 1 s in several transition 
region line including the intense Si IV line at 1394 A. 

Such high cadence spectral data has the potential to provide diagnostic in- 
formation on the plasma's electron temperature and density. To do so requires 
transient ionization considerations. As a test on the diagnostic potential we 
analyze one of the flares observed by SMM, making use of the Atomic Data and 
Analysis Structure (ADAS; Summers, 200E@) framework, which is a collection 
of fundamental and derived atomic data, and codes that manipulate them.. 

In Section 2, we briefly outline past observations and findings on the temporal 
and spatial correlation between HXR and UV emissions during solar flares and 
the high cadence UVSP data used in this paper, while in Section 3, we outline 
the atomic background, give details on a full time-dependent calculation and 
comparing it to the results obtained for the radiative power assuming steady 
state in Section 4. Here, we only consider an atomic model, i.e. we assume that 
the plasma is instantaneously heated from a cold plasma at 2 — 3 x 10 4 K to 
temperatures 3 — 5 x 10 5 K. We do not consider movement of plasma nor a model 



1 http: / /iris. lmsal.com/ 
2 http:/ /www. adas.ac.uk 
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atmosphere, as was done for comparison with jet-like features by Roussev et al. 
(2001), Doyle et al. (2002). 

Most radiated power loss calculations are done on the basis of a steady-state 
ionization balance assumption, which is relatively simple and is often used in 
circumstances which are not fully justified. Non-equilibrium ionization occurs 
when the plasma electron temperature or electron density changes on a timescale, 
Tpiasma, shorter than the atomic ionization stage fractional abundance relaxation 
timescale, r z (where z is the ionization stage of a specific element), or when the 
time for plasma transport across a temperature or density scale length is shorter 
than T g . In the present work, the first case is considered. Following McWhirtcr 
(1965), for an ionizing plasma, r 2 represents the relaxation time for a plasma 
to reach ionization equilibrium. Transient ionization was considered by Doschck 
and Tanaka (1987) for 20 x 10 6 K flare phase while Raymond (1990) looked at 
its effect for a micro-flare heated corona. In each instance, these authors noted a 
flux enhancement resulting from transient ionization considerations. Due to the 
general lack of high time resolution UV line flux observations, whether from flares 
or active regions, this aspect is not generally appreciated. In the present work 
the importance of non-equilibrium ionization is shown comparing the O v 1371 
A line emission enhancement (with respect to equilibrium) predicted by a model 
with the observation of a flare which occurred on 2 November 1980 (Poland et 
al. 1984). 

2. Temporal and Spatial Correlation Between HXR and UV 
Emissions During Solar Flares 

2.1. Past Observations 

The relationship between HXR and UV emissions during flares has been inves- 
tigated since the early 70s (Kane and Donnelly (1971) and references therein), 
using broadband EUV emission. Uatcr Poland et al. (1984) and Cheng et al. 
(1981, 1988) found a strong temporal correlation between HXR emission at 25- 
100 kcV and the O v 1371 A emission, which peak at 25 cV (« 2.9 x 10 5 K), 
during the impulsive phase of solar flares. They used simultaneous observations 
obtained by the Hard X-Ray Burst Spectrometer (HXRBS) and Ultraviolet Spec- 
trometer Polarimetcr (UVSP) on the Solar Maximum Mission (SMM) satellite. 
The HXRBS had a temporal resolution of 0.128 s, but no spatial resolution, 
so it gave HXR emission from the entire Sun, while the O v emission was 
from the pixels showing flare enhancements. Poland et al. (1984) studied the 
temporal correlation between the HXR and O v emissions using a sample of 
17 flares, finding a well defined relation between the two emissions for a given 
flare. However, this relation varies from one flare to the next. They analysed 
the observations in the context of a single loop model heated by a high-energy 
electron beam. The beam is accelerated down to the foot-points in the chromo- 
sphere. Here the beam collides with the cool chromospheric plasma and yields 
HXR brcmsstrahlung emission, which heats the plasma through the Coulomb 
collision with the electrons. The plasma heated to transition region temperature 
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Figure 1. The UVSP and HXRBS data for a flare which occurred on 2 November 1980 (see 
Poland et al. 1984). The upper plot shows the O V 1371 A light curve with the individual pulses 
labeled 1 to 10. The solid line represents the observation, the dotted line is the fit performed 
using the simple function given by Equation [6] and the red line is the sum of all the bursts 
(see Section 4). The lower plot displays the O V line light curve together with the hard X-ray 
emission to show the correlation between the O V bursts and the hard X-ray pulses. 
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gives increased O v emission. They concluded that there is a strong physical 
relationship between the processes which give rise to the non-thermal HXR 
emission and the apparently thermal UV emission. However, while they could 
predict with their models the general behaviour of O V emission compared to the 
HXR, they did not reproduce the magnitude and detailed behaviour. A close time 
correspondence between the HXR and O v emissions was found also by Cheng et 
al. (1988), who analysed five flare events. Their cross-correlation analysis showed 
that the fast features of O v generally lag behind the corresponding HXR features 
by 0.3-0.7 s. A revisit of the relationship between the UV and HXR emissions in 
flares was done by Warren and Warshall (2001). They studied nine flares, using 
data from it TRACE, it Yohkoh Hard X-Ray Telescope and Burst and Transient 
Source Experiment (BATSE) on the Compton Gamma-Ray Observatory, with 
a cadence of 2-3 s. For five of the nine flares their observations allowed the 
spatial analysis of the HXR sources (with a spatial resolution of 8"), finding for 
2 flares a spatial separation between the HXR and UV sources. Further detailed 
studies on both temporal and spatial correlation between the HXR and UV 
emission in flares have been done by Alexander and Coyner (2006) and Coyner 
and Alexander (2009). They incorporated the improved spatial and temporal 
resolution of it TRACE (focusing on the 1600 A observations) and RHESSI. 
The latter allows the examination of spatial HXR structures with a resolution 
of 2", while it permits a time resolution of tens of milliseconds. However, the 
temporal analysis was limited by the it TRACE cadence of 2 s. They investigated 
the spatially integrated UV and HXR time profiles. In addition, they compared 
the individual UV sources (or group sources) with both the overall HXR time 
evolution and the individual impulsive X-ray bursts. Finally, they performed a 
spatial association between the spatial resolved UV and HXR sources. Alexander 
and Coyner (2006) focused on a particular flare, which occurred on 2002 July 
16 and found a good temporal correlation between the spatially integrated UV 
and HXR emissions (as found by previous authors), but different temporal and 
spatial correlation for the individual sources. Three main conclusions arose from 
their work: 

(i) the strong overall correlation between the UV and HXR emissions suggests 
that these two different emissions must result from the same energy release 
process or a strictly related processes; 

(ii) the separation of the UV emission into individual sources provided examina- 
tion of the various contributing components to the overall temporal behaviour 
and suggests a more complicated overall picture of solar flares than the single 
loop models of previous work; 

(iii) the lack of spatial correlation for some UV and HXR sources suggests a 
three-dimensional magnetic topology. 

Coyner and Alexander (2009) improved the analysis of Alexander and Coyner 
(2006) investigating the behaviour of temporal and spatial UV and HXR emis- 
sions in a large sample of solar flares. They found that the HXR emission 
characterized by multiple impulsive bursts, that each individual burst corre- 
sponds to independent spatially disconnected sets of HXR foot-points. These 
sources are localized and are co-spatial with some of the UV sources. However, 
while the HXR sources are compact and localized, the UV emission is more 
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diffuse and persists throughout the entire flaring region. They noticed that 
some UV sources were not co-spatial but were co-temporal with respect to 
the HXR sources and some other UV sources showed neither co-spatial nor 
co-temporal emission compared to HXR sources. They explained the temporal 
but non-spatial correlation between some UV and the HXR sources, suggesting 
a three-dimensional magnetic field morphology, again more complex than single- 
loop models. Concerning the spatially and temporally uncorrclated UV emission, 
they pointed out that these signatures indicate two distinct types of UV response, 
one consistent with the HXR emission and one of a likely thermal origin (which 
may be indicated by the significant correlation with soft X-ray emission at 6-25 
keV instead of the HXR at 25-100 keV). They concluded that in solar flares the 
structure of the magnetic field must be a complex three-dimensional structure 
and that there arc multiple processes involved in the UV emission production. In 
summary individual bursts in the HXR time profile are related to disconnected 
sets of HXR foot-points and that a probable origin of the HXR emission is related 
to electrons (and ions) accelerated in the corona, traveling to the foot-points in 
the chromosphere. Here these electrons heat the plasma up to transition region 
temperature producing the UV emission. 

There are three types of relationship between UV and HXR emission: 

(1) co-temporal and co-spatial UV and HXR emissions; 

(2) co-temporal but non co-spatial UV and HXR emissions; 

(3) neither co-temporal nor co-spatial UV and HXR emissions. 

The first type of relationship, that is the existence of localized HXR emission 
co-spatial with temporally correlated UV emission sources, is consistent with 
models where both emissions result from different portions of an accelerated 
electron population. The second type of relationship, that is spatially separated 
but temporally correlated UV and HXR emissions, is consistent with a picture 
of a spatially complex magnetic geometry for solar flares, different from a single 
loop model. The third type of relationship, that is spatially and temporally 
uncorrelated UV and HXR emissions, is consistent with the presence of two 
distinct mechanisms each producing UV emission (e.g. one consistent with HXR 
origin and dependent on it and one of a likely thermal origin and independent 
of HXR emission). 

2.2. The Data Used Here 

The Ultraviolet Spectrometer and Polarimeter (UVSP) is described by Woodgate 
et al. (1980), while details on the Hard X-ray Burst Spectrometer (HXRBS) is 
given by Orwing et al. (1980). For the flare observations discussed here, the 
entrance aperture was 10", with an exit slit of 0.3 A centered at O V 1371 A. 
In response to an alert from another onboard instrument, the UVSP began 
making rasters covering 30" x 30". In some of the earlier UVSP/HXRBS flare 
observations, the UVSP data was taken with a 1.3 s cadence and in order to 
compare the two datasets, the UVSP data was interpolated at the finer HXRBS 
time resolution of 0.128 s. However, in most of the later observations, both 
datasets were observed using a 0.128 s cadence. In order to improve the S/N, 
several of the bright UVSP pixels were added. 
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Of the various papers which reported the almost simultaneous UV/hard X- 
ray bursts, Cheng et al. (1988) did a detailed filtering with a cross correlation 
analysis. These authors established that the O V 1371 A line lagged the hard 
X-ray burst by 0.3 to 0.7 s, and was flare dependent. Here, we use this data 
to show the diagnostic potential of high cadence UV spectral data. In Fig. 1, 
we reproduce data from Poland et al. (1984) of a flare which occurred on 2 
November 1980. 

3. Atomic Background 

To predict the line emissivities and ion population densities, the generalized 
collisional-radiative (GCR) theory as implemented in ADAS, the Atomic Data 
and Analysis Structure (McWhirter and Summers, 1984; Summers, 2009), is used 
in this work. Each ion in an optically thin plasma is described by a complete 
set of levels with collisional and radiative couplings between them. All radiative 
and electron collisional processes are included in the model, except for photon- 
induced processes. State-resolved direct ionization and recombination to and 
from the next ionization stage arc also taken into account. In addition, due to 
the much shorter relaxation timescales for the excited states (not metastables) 
with respect to the ground and metastable states, it is assumed that the dom- 
inant populations lie in the ground and metastable states of the ion. Over the 
considered temperature and density regime, the ratio in the population of the 
metastable state 2s2p 3 P compared to the ground state 2s 2 1 S of O v is typically 
less than 10 -2 . Hence the dominant population is concentrated in the ground 
state. 

In a time-dependent plasma model the line emissivity is no longer a unique 
function of the local temperature and density conditions but depends on the 
past history of the temperature, density and state of ionization of the plasma. 
Therefore, the assumption of ionization equilibrium in calculating the ionization 
balance is not appropriate and time-dependent fractional abundances must be 
determined. The time dependence of ionization stage populations, N^- z \ leads to 
the following equation: 



where the presence of metastable states is neglected because of the previous 
consideration. S and a are the collisional-dielectronic ionization and recom- 
bination coefficients. They give the contribution to the growth rates for the 
ground state population, due to the effective ionization, which includes direct 
and excitation/auto-ionization contributions, and the effective recombination, 
which includes radiative, dielectronic and three-body contributions. The values 
of these coefficients, currently within the ADAS database, have been obtained 
following the GCR approach as described in Summers et al. (2006). The solution 
of Equation [1] is such that the number density of the element of a nuclear charge 
Z, is equal to X^z=o ■ ^ ne time-dependent fractional abundances 

are calculated as follows: the code derives the solution for a range 



dt 



JV e [S ,(z_1) iV (2 " 1) + (5* (z) + a (z) )A^ (z) + a (z+1) A^ (z+1) ] 



(1) 
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of fixed plasma electron temperature and density pairs, starting from an initial 
population distribution N^ z '(t = 0)/N^ (see Section [4]) using an eigenvalue 
approach. 

The line cmissivity is given by: 

£j ^=A cI ^7V c 2 Gg 4 (T c ,iV c ,i) (2) 

where A e \ = /TVh is the abundance of the element Z relative to hydrogen, 
N R /N C is tabulated by McWhirter et al. (1975) and Gf\ t (T c , 7V c t ) i s the time- 
dependent contribution function defined as: 

G%(T c ,N Cl t) = VECf^^^+VEC^^^ (3) 

where 

A^F^^ (4) 

Aj^i^ 16 ^ (5) 

where Aj-^ is the transition probability and 7V( Z+1 ) are the population 
densities of the ground states of the ions of charge z and z + 1 and J r ( cxc ^ z ) and 
are the effective contributions to the population of the upper excited 
state i. Further details on the GCR may be found in Lanza et al. (2001) 

4. Burst Model and Conclusions 

The main goal of this work is the analysis of the effects of a plasma not being 
in ionization equilibrium due to rapid heating that is faster than the ioniza- 
tion/rccombination timcscalcs. The issue is addressed considering each burst (10 
in this case, see Fig. 1) as an independent entity, which is assumed to originate 
as a consequence of a plasma heating at the foot of a loop in a multiple loop 
structure. Whether this energy release results from a non-thermal electron beam, 
and/or a conduction front is not the subject of this paper. The calculations 
performed for this work consider only the plasma response to sudden heating 
with no mass motions allowed. 

Using the above considerations on the strong physical link between the HXR 
emission and the O v 1371 A emission, the initial conditions for the plasma 
evolution are stated (see Table 1). For each burst it has been assumed that at 
time t = 0, the plasma is in ionization equilibrium at an electron temperature 
in the range 3 — 5 x 10 4 K and the electron density in the range 5 x 10 10 — 
5 x 10 11 cm -3 (typical values of the solar chromosphere). Using these initial 
conditions, the whole population of oxygen lies in the first four ionization stages 
(O I - O iv) and is concentrated essentially in the second and third ionization 
stages (O II and O in), as shown in Table[TJ This Table lists only the range of the 
initial temperatures and densities considered because these initial conditions do 
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not significantly affect the late plasma evolution when the spectral lines observed 
are emitted. They depend strongly only on the final temperatures and densities. 
Then, to simulate a flare burst, the temperature of the plasma is instantaneously 

Table 1. Initial conditions for the plasma 
evolution. At t = the fractional abundances 
have been calculated in equilibrium condi- 
tions, assuming T c = 3 — 5 X 10 4 K and 
N e = 5 X 10 10 - 5 X 10 11 cm- 3 . The ta- 
ble shows the fractional abundances for the 
most populated ionization stages in the stated 
(T e ,N e ) conditions. 





N c \ T e 


3 x 10 4 K 


5 x 10 4 K 


5 x 10 10 cm" 3 


O 11=0.98 
O iii=0.01 


O 11=0.28 
O iii=0.71 


5 x 10 11 cm -3 


O 11=0.98 
O iii=0.01 


O 11=0.25 
O iii=0.74 





raised by a factor w 10 (from around a few times 10 4 K to around a few times 
10 5 K) with the bursts of Fig. 1 reconstructed following several steps, as listed 
below: 

(i) Each observed burst has been fitted using the function: 

f(t) = N oh B R e -(^) (6) 

where B = (t — ts)e 1 Rr 1 / r oh . r ob is an exponential time decay, t$ is the start 
time of the burst, N ob is the peak line enhancement value of the burst and R 
controls the width of the burst. This analytic expression has been used in the 
past in the context of the "smooth-burst model" (Doyle et at 1991; Kellett and 
Tsikoudi (1999) and references therein) to estimate flare intensity profiles for 
different types of stars. Here this model has been used to reconstruct the profiles 
of the observed solar bursts and estimate their decay times (r ob ) and peak values 
(iV ob ). Table [2] lists the values derived fr om this fitting procedure for each burst. 

(ii) The modelled relaxation timescales, r th of the ionization/recombination 
processes for O v have been calculated for different pairs of electron temperature 
and electron density using the relation given below (McWhirter, 1965). This 
relaxation time characterizes the time to reach ionization equilibrium and is 
defined as the reciprocal of the product of electron density and the sum of 
ionization and recombination coefficients relative to the considered ionization 
stages: 



Ne ( S ( z ^ z+ i) + a (z+l->-z)) 
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Table 2. Values of de- 
rived decay time (T ob ) and 
peak line enhancement fac- 
tor (N° b ) for the ten ob- 
served bursts derived via 
fitting the function given 
in Section 4. 



Burst 


r ob (s) 


N° b 


1 


19.41 


2.46 


2 


7.05 


2.76 


3 


6.41 


3.11 


1 


6.22 


3.25 


5 


4.20 


2.94 


6 


4.31 


3.33 


7 


4.05 


3.68 


8 


8.38 


2.65 


9 


8.03 


2.36 


10 


20.39 


2.19 



where z = 4 is the ionization stage with which the O v 1371 A line is concerned, 
^(z-fz+i) Quotes the effective ionization coefficient and a( z+1 ^ z ) the effective 
recombination coefficient. 

(iii) A theoretical peak line enhancement factor value 7V th represents the 
theoretical line enhancement at the peak of the burst with respect to the equi- 
librium value. An example is shown in Fig. [5J This figure plots the transient 
line contribution function G*.^_ >i .(T e , N e , t) as a function of time compared with 
the equilibrium line contribution function G^ fc (T e , N e ) for two pairs of electron 
temperature and electron density values, which give the same r th but different 
line enhancement factor A th . 

(iv) The theoretical relaxation timescales r th and peak line enhancement 
factor values N th have been derived for a set of (T e ,N e ) pairs for the calculation 
of the transient fractional abundances and equilibrium fractional abundances, 
respectively. Over-plotting contours with the observed values (r ob and iV ob ) 
listed in Tabic [H the final electron temperatures and electron densities for the 
time evolution of each burst are derived. Figure [3] displays an example for bursts 
4 and 7. The lines with r constant (solid lines), corresponding to the decay 
times of the fourth and seventh bursts, intersect the lines with N constant 
(dashed lines), corresponding to the peak values of the two considered bursts. 
The intersection points give the electron temperatures and electron densities for 
the time evolution of these two bursts. The same procedure has been applied to 
the other bursts. The results are summarized in Figure 01 which shows a contour 
plot of the theoretical relaxation timescales as a function of electron temperature 
and electron density. 

The over-plotted black lines are the contours which correspond to the observed 
decay times of Table [3] The white points represent the intersection with the 
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Figure 2. Contribution functions of O V 1371 A as a function of time for reconstructing 
the second burst. The two solid straight lines are the equilibrium contribution functions 
calculated for two pairs of electron temperature and electron density values (the blue one 
is G -Z+kiTe,!, ^Ve,l) while the green one is G°^ fc (T e ,2, N c ,2) )■ The two curves which change 
with time represent the transient contribution functions calculated for the same two pairs of 
electron temperature and electron density. They are characterized by the same relaxation time 
r th _ r ob _ 7 Q5 s but diff erent pea k values, 7V th (T Ci i, AT ,i) = 2.76 (which is equal to the 
observed value) and JV th (T 8 ,2, iV e , 2 ) = 1.24. 



iV=constant lines for each burst. The whole set of points lies in a region of 
quasi-constant electron temperature (T w 2.6 x 10 5 K), while there is a larger 
spread in the electron density domain (N c w 2.6 x 10 9 — 1.1 x 10 10 cm -3 ), see also 
Figure [S] This figure shows that near to the highest peak (i.e. burst 7, Figure 1), 
the density is high while the relaxation time has one of the lowest values. The 
longest relaxation times are related to the first burst (burst 1) and last burst 
(burst 10). The latter is evidence of the slow decrease in the intensity of O V 
after the flare event due to the plasma relaxation until it reaches equilibrium 
conditions. 

Figure [5] shows the contribution functions as a function of electron tempera- 
ture. Note the much larger contribution function resulting solely from transient 
ionization. The dashed lines represent the instantaneous contribution functions 
calculated at the time where each burst peaks using the respective electron 
densities and the solid lines are the equilibrium contribution functions calculated 
at the same densities. The non-equilibrium contribution functions peak between 
2.6 — 2.9 x 10 5 K while the equilibrium contribution functions peak at 2.1 x 10 5 K. 
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Figure 3. Lines with r constant (solid lines), corresponding to the decay times of the fourth 
and seventh bursts. The dashed lines have N constant, corresponding to the peak values of 
the two considered bursts. The intersection points give the derived T e and We 



The final model to reconstruct the time evolution of each burst has been built 
up considering the following relation between the line emissivity calculated in 
non-equilibrium conditions at a fixed pair of electron temperature, and electron 
density and equilibrium line emissivity calculated for the same (T e ,N e ) pair: 

e^ k (T^,N oA ,t) _ G^ k (T 0tU N 0tU t) , . 

e^ k (T e , u N e;i ) G^ fc (T e)i ,JV eii ) * ' ' W 

where i is the burst number. Equation [8] provides the enhancement factor due to 
transient ionization with respect to the equilibrium conditions, see the red line 
in Figure [T] which is the sum of all the bursts. 

In the previous section we used initial temperatures and densities typical 
of the solar chromosphere, in order to check the sensitivity of these values we 
performed additional calculations for one burst, in this instance burst 3. This 
burst showed an observed enhancement factor of about 3 and a decay time of 
about 6 s. It is possible to reproduce these values with a final JV e « 7 X 10 9 
cm~ 3 and T c w 2.6 x 10 5 K. Increasing the final temperature to 1 x 10 6 K (with 
the initial temperature still in the range 3 — 5 x 10 4 K), and a final electron 
density in the range 10 9 to 10 13 , the enhancement factors becomes 1 x 10 5 to 
1 x 10 6 instead of 3, while the decay times decreases from about 2 s to less than 
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9.0 9.2 9.4 9.6 9.8 10.0 10.2 10.4 
log(density) (cm -3 ) 



Figure 4. A contour plot of r th calculated using Equation 171 The black lines are the constant 
rs (=T ob ) obtained by the fitting procedure. The white stars are the derived values of T e and 
N e for each burst, determined as in Fig. [3] The lower contour plot is an enlargement of the 
upper black rectangle showing the position of each burst labeled by numbers from 1 to 10. 
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Figure 5. The derived electron temperature, electron density and relaxation time as a function 
of the time at which each burst has its peak. 
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Figure 6. O V 1371 A contribution functions as a function of electron temperature. The 
dashed lines are the instantaneous contribution functions at the time where each UV burst 
peaks. The corresponding electron density is given for each burst. The solid lines represent the 
contribution functions calculated using equilibrium ionization balance at the same electron 
density. 



0.01 s as the density increases. Changing the initial temperature to 1 x 10 5 K, 
only results in the final values changing by a few percentage. This means that 
in order to re-produce the observed enhancement and decay time, these initial 
conditions are appropriate. 

The conclusion from the results given above is that even a small departure 
from ionization equilibrium can result in an enhancement of the O v 1371 A line 
emission compared to the intensity calculated in ionization equilibrium. The 
present model reproduces the time evolution of each burst, although some issues 
need further investigation, e.g. it is noted that all the derived points in Fig. 1 
lie in a very small temperature range. Thus more observations, together with a 
statistical analysis are needed to investigate whether this temperature behavior 
is the same for other flares. 

This study assumes a Maxwellian electron distribution, which is justified by 
the substantially shorter average electron-electron energy relaxation time com- 
pared with the observed plasma development and with the theoretical ionisation 
stage (ground and metastables) relaxation times. However, the differential vari- 
ation of electron-electron relaxation, especially to higher energies could allow a 
substantial cohort of fast electrons to confuse interpretation. It is unlikely that 
a sufficient set of line observations will be available for the solar atmosphere 
to explore this, however, it is hoped that parallel studies in tokamak divertor 
transients will clarify experimentally the non-Maxwellian spectral indicators. 

What this work clearly shows is that observed transition region line in- 
tensity increases on a sub-second time-scale can be explained via transient 
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ionization. The response time is electron density dependent, explaining the re- 
ported O v/HXR time delays of Cheng et al. (1988). This work has implications 
not only for high cadence flare observations, but high cadence observations of 
other dynamic solar features (Doyle et al. 2006, Madjarska et al. 2009). There 
are currently no spectrographs capable of sufficiently high cadence observations, 
although the selected IRIS mission (launch date 2012) will be capable of obtain- 
ing such data for a range of transition region lines including Si iv 1394 A and 
O IV 1401 A. The contribution functions for these two lines behave differently 
theoretically under transient ionization, with Si iv 1398A more like O v 137lA. 
Both of these lines are in the ground state spin system. The O IV 1401 A line's up- 
per state is in the quartet (mctastablc) spin system. Transient ionisation under- 
fills the population structure of the alternate spin systems from the ground. By 
contrast the ground spin system population structure is enhanced with the expo- 
nential factors in the rate coefficients largely unwound because of the increased 
electron temperature. 
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